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Abstract 

We study the de-Haas van Alphen oscillations in the magnetization of 
the Hofstadter model. Near a split band the magnetization is a rapidly 
oscillating function of the Fermi energy with lip shaped envelopes. For 
generic magnetic fields this structure appears on all scales and provides 
a thermodynamic fingerprint of the fractal properties of the model. The 
analysis applies equally well to the two dual interpretations of the Hofs- 
tadter model and the nature of the duality transformation is elucidated. 

The Hofstadter model [1] describes non-interacting (spinless) electrons mov- 
ing in the plane under the combined action of magnetic field and a periodic 
potential. The model is a paradigm for quantum systems with singular con- 
tinuous spectra [2, 3]; Chern numbers in the integer quantum Hall effect [4]; 
Fractal quantum phase diagrams [5]; Quantum integrable models [6] and Non- 
commutative geometry [7]. It has been realized experimentally in [8] and mea- 
surements of the Hall conductance phase diagram are in good agreement with 
theoretical predictions of [4]. 

Previous studies of the thermodynamics of the Hofstadter model have fo- 
cused on minima of the ground state energy [9]. Our aim is to examine the 
possibility that the de Haas- van Alphen oscillations, which are a basic tool in 
the study of the geometry of Fermi surfaces of metals [10], also provide a tool 
for studying its non-commutative analog. We show that when the magnetic 
flux per unit cell is close to rational multiple of the flux quantum, the magne- 
tization oscillates as a function of the chemical potential in the vicinity of the 
bands of the rational flux, as in the de Hass-van Alphen effect. As the deviation 
from rational flux becomes smaller, the frequency of the oscillations increases, 
so that in the limit of the deviation tending to zero the oscillations fill an area 
bounded by a limiting envelope. The envelopes have the following universal 
features: They are smooth functions of the chemical potential except for loga- 
rithmic pinching at a single point, associated with a logarithmic divergence of 
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the density of states, and they are pinched linearly at the edges of the bands. 
These properties give the envelopes suggestive lip-like shapes. 

It follows from these results that at magnetic flux ratios which are irrational 
numbers well approximated by infinitely many rationals numbers, magnetic os- 
cillations bounded by lip-shaped envelopes occur on all energy scales. In this 
sense, the magnetic oscillations provide a fingerprint of the fractal spectra of 
the model. 

A second issue of the Hofstadter model that we discuss is duality. The 
Hofstadter model approximates the Schrodinger equation in two dual limits: 
When the magnetic field is a weak perturbation of a Bloch band and remarkably 
also in the dual limit when the magnetic field is so strong that the periodic 
potential is a perturbation of the lowest Landau level. As we shall see, the 
thermodynamic properties of the dual limits are related by an interesting duality 
transformation of the thermodynamic potential, Eq. (3) below. We examine the 
consequences of the duality for the magnetization and the Hall conductance in 
some detail. 

The spectral analysis of the Hofstadter model reduces to studying a family 
of one dimensional problems parameterized by a conserved quasi- momentum £. 
For the sake of concreteness and simplicity we shall consider the case where 
the periodic potential has square symmetry. In this case the corresponding 
Hamiltonian is: 

= 1>(j - + 1) + 2 008(2^ + 0^0'). (!) 

It is known that H describes the splitting of a Bloch band by a weak transverse 
magnetic field, and also the splitting of a Landau level by a weak periodic 
potential, <p> and tp^ 1 are the values of flux through a unit cell, measured in 
units of quantum flux $o = hc/e, in the split Bloch and split Landau cases, 
respectively. Note that ip is a dimensionless quantity. 

Although H(ip, £) depends analytically on ip its spectral properties are sensi- 
tive to the flux: When ip is a rational number the spectrum is a finite collection 
of energy bands, while for (almost all) irrational fluxes it is a Cantor set of zero 
Lebesgue measure [3]. 

Let fi(,(T, n,<p) be the thermodynamic potential per unit cell of a split Bloch 
band. By gauge invariance, time-reversal and electron-hole symmetry, Q b satis- 
fies 

Q b (T,ij,,ip) = Qb(T, n, —ip) 

= il b (T,fj,,<p+l) =fj, + il b (T, -n,<p) (2) 

It follows that the thermodynamic properties for all tp and \i are determined by 
those for < tp < 1/2, [i < 0. yu = corresponds to the half-filling point of 
electron-hole symmetry. The phase diagram of the split Bloch band is displayed 
on the left hand side of Fig. 1. The phases are labelled by the values of the 
Hall conductance, whose relation to f2 is given in Eq. (4) below. The phase 
diagram displays the symmetry properties shown in Eq. (2). The phase diagram 
is periodic in (p, sec Eq. (2); a single period is given in Fig. 1. 
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Figure 1: The quantum phase diagrams of a split Bloch band (left) and split 
Landau level (right). The horizontal axis is the chemical potential and the 
vertical axis is < ip < 1 on the left and < (p^ 1 < 2 on the right. The colors 
represent the quantized Hall conductances. The scheme is such that cold colors 
correspond to negative integers and warm colors positive integers. Zero is white. 
The red on the right hand side of the right diagram corresponds to unit Hall 
conductance of a full Landau band. 

It is straightforward to arrive from the the thermodynamics of the split Bloch 
band to that of the split Landau level by observing that the energy spectra of 
the two models are the same. The difference is that the number of states per 
unit area in a split Bloch band is a constant independent of the magnetic field, 
whereas the number of states per unit area of a Landau level is proportional to 
the magnetic field. Therefore the thermodynamic potentials of a split Landau 
level, fli, and split Bloch band are related by 

tti(T,H,p) = ipttb{T,n,ip~ l ) (3) 

This is a duality transformation: It is symmetric under the interchange b < — ► I. 
It implies that the thermodynamics of the split Bloch band determine the ther- 
modynamics of a split Landau level and vice versa. It follows from Eq. (3) that 
fli is not periodic, although fif, is, simply because of the factor p on the right. 
A part of the phase diagram of the split Landau level is shown on the right hand 
side of Fig. 1, and it is evidently aperiodic. 

To derive the relation between the magnetization and Hall conductances of 
the two models recall that the filling fraction p, the magnetization per unit area, 
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in, and Hall conductance [11], a, are given by 

fdfl\ 1 fdn\ e 2 (dp 



U*V "— *o V<W„' " h v^;/ (4) 

The magnetization and the Hall conductances of the two models are therefore 
related by: 

mii^T,^- 1 ) = -^-n b (fj,,T,<p)-<pm b (fJ,,T,ip); 
e 2 

at(H,T,<p ) = — p b (p,T,tp) - <pa b (p,T,(p) (5) 

The duality relation of the Hall conductances is a generalization of a relation 
that follows from the Diophantine equation in [4]. 

The magnetization is related to the Hall conductance by the Maxwell relation 



O y = io(!X=^- (6) 



It follows that whenever the Hall conductance is quantized, the magnetization 
has a quantized slope that is an integer multiple of l/$o- In the Hofstadter 
model a b and 07 are quantized in the gaps [4, 12]. It follows that the magneti- 
zation has quantized (and therefore also constant) slopes in the gaps. 

We now turn to a more detailed study of the magnetization in the gaps at 
T = 0, starting with m b . For a fixed rational ip this is a finite collections of 
lines. From Eqs. (2), (4) it follows that m b is a symmetric function of p. 

When p = p/q, with p and q relatively prime, the magnetization can be 
computed numerically. In this case the potential in Eq. (1) is periodic with 
period q and the spectral analysis of H(ip,£) further reduces to the study of 
q x q matrices H(<p, £,77) labelled by a second Bloch (quasi) momentum 77. The 
magnetization m b (p,tp) in the gap of a split Bloch band can be shown to be 
given by 

(7) 

e are the eigenvalues and u the corresponding eigenvectors of the matrix. > , < 
refer to states above and below the gap. Our numerical results are based on 
this equation. 

Our investigation of the magnetization in the gaps bears an interesting rela- 
tion to the work of [9] , who study of the energy E per unit cell of the Hofstadter 
model. E is the Legendre transform of fl, 

E(p,p) = pp-fl(p,tp) , (8) 

and [9] focused especially on its minima for fixed p. The properties of the 
Legendre transform imply that maxima of as a function of <p, i.e., points where 
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m = and a > 0, are (local) minima of E. It follows from our semi-classical 
analysis (see below) that E(p, tp) has a minimum in every gap for ip — 1/q, see 
for example Fig. 2. It also follows from Eq. (8) that E has a cusp as a function 
of tp at the gaps, with tangent slopes given by the limiting values of m at the 
ends of the gap. 

-4> m b 




Figure 2: The graph of the magnetization m& of a split Bloch band with tp = 
1/101. When tp tends to zero the magnetization has a limiting envelope denoted 
byL fc ( M ,0). 

The Hofstadter model can be viewed as the result of Peierls substitution [13] , 
which stipulates that in a weak magnetic field the classical Hamiltonian 

H B {Z,v) =2cos£ + 2cos?? (9) 

is quantized by imposing [£, 77] = 2tt i tp. Fluxes near tp = can then be analyzed 
by semiclassical methods where the role of H is played by 2irtp. The classical £ — r\ 
phase space is the two dimensional torus. Since its area is (27r) 2 , quantization 
of the torus is consistent provided the number of states, l/tp, is an integer. For 
this reason we shall assume that tp = 1/q. 

The classical Hamiltonian Hb{S,,t]) is a Morse function with one minimum, 
(at energy —4), one maximum, (at energy 4), and two saddle points (at energy 
zero). It follows that the classical trajectories with E ^ are contractible and 
bound a disc like domain. Only the level set E = winds around the torus. 

Let 

S(E) = [ dl; A drj (10) 

JH B <E 

be the classical action associated with a trajectory of energy E. By the Bohr- 
Sommerfeld quantization rule, the negative part of the spectrum is given, to 
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leading order in ip, by those energies E n for which 

S(E n ) = (n + -y)<p, 7 = 1/2, n=l,...,L«/2j (11) 

This gives a good approximation to the spectrum provided that ip is small and 
E is far from the separatrix. More precisely, one needs tp <C \E\. The energy 
E n of Eq. (11) is then an approximation to the n-th energy band which is 
exponentially localized in energy near E n . Near the separatrix one finds wider 
bands which we do not study [14]. We therefore assume below that /i <C —(p. 
The magnetization for /i>^ then follows by Eq. (2) . The region which we do 
not study, \/j,\ — 0{<p), becomes arbitrarily small in the limit p — > 0. This is a 
limit we eventually take. 

Fixing the chemical potential in the n-th gap with fj, < 0, the semiclassical 
thermodynamic potential at T = is 

n— 1 n— 1 

n b = HP~Yl E i =Vn<p-tpJ2 S ^ (0' + (12) 
j=o j=0 

By the second Euler-Maclaurin formula, the sum can be approximated by an 
integral (to 0(1 /q 2 )), 

flb(p,p) = finip — / S~ 1 (x)dx, (i < 0. (13) 
Jo 

The magnetization in the n-th gap is therefore 

m b = -r-j- (a*- S-^ntp^n (14) 

and vanishes when /j, = S l_1 (n< ) 5). Since 7=1/2 this happens in the middle of 
the gap, a property also evident in Fig. 2. 

Since the slope of nib in the n-th gap is n/$o, the magnetization increases in 
this interval from —n(E n - £„_i)/2$ to +n(E n - £„_i)/2$ - The extremal 
value Lb of the magnetization in the n-th gap is therefore 

L b = ±^-S- 1 '(n V ). (15) 

Fixing the chemical potential /x and letting <p — > 0, the product n<p approaches 
the limiting value 5(/z), see (11). The envelope approaches 

^-"'^s^b- (16 » 

The graph of L b (fj,, 0) is shown in Fig. 2 overlayed with numerical results for the 
magnetization that follow from the (exact) Eq. (7). This graph is an analog of 
the de Haas- van Alphen oscillations in metals [10]. 

L b (n, 0) is the mother of all lips: We shall argue below that the lips at other 
(rational) fluxes are deformations of it, in both the the split-Bloch and split- 
Landau versions of the model. The lips associated with the split Landau level 
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£z(/z, 1/n), with integer n, are particularly simple deformations of the mother 
of all lips. This can be seen as follows: From Eq. (5) 

Li(fi, 1/n) = — n b (fi,oo) -nL b (n,0) (17) 
'I'o 

Since £l b is a convex, and therefore continuous, function of fi, Eq. (17) describes 
a scaling of L b ([i,0) and a shifting of its baseline, see Fig. 3. The special case 
rni(/x, oo ) (n = 0) is different: The lip is scaled to zero width, and m\(ji, oo) is 
a continuous function of \i. 
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Figure 3: The central curve gives the limiting magnetization, m;(/i, ip = oo) and 
the upper and lower curves form the lip £;(/x, ip = 1) 

Although the actual shape of the lip L b (/j,, 0) cannot be expressed in terms of 
elementary functions, its qualitative features can be readily understood and are 
determined by the critical points of the band function Hb{^, ??)■ Let 5(x denote 
the distance of \x from a band edge. Since Hb is non-degenerate, the action, 
S(/i) is linear in <5/x near the maximum and the minimum — as for the harmonic 
oscillator. This says that the lip terminates linearly: 

S(m)~|<H L b ((M,0)^±5^. (18) 

Near the separatrix at E — 0, which passes through the two saddle points of 
H(£,r]), the density of states diverges logarithmically and L b vanishes logarith- 
mically 

S( M )-S(0)~ M log| M |, L b (n 1 0)~±(\og6ti)- 1 . (19) 

Consider now the magnetization for tp close to a rational p/q. Since the gap 
edges are Holder continuous of order 1/2 in ip [15], the q — 1 gaps (q — 2 for 
even q) persist for nearby fluxes and so does the magnetization in these gaps. 
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Figure 4: Plot of the magnetization as function of the chemical potential for 
flux that is 40/121 of the unit of quantum flux through the unit cell of the tight 
binding Hofstadtcr model. The magnetization is plotted for values of chemical 
potential that correspond to a spectral gap. 

The q bands, however, splinter. This breakup can be described by semiclassical 
methods as follows [2]: For rational flux <p — p/q, the q x q matrix H(tp, £,77) 
gives rise to q magnetic Bloch bands Ej(£,rf), j = l,...,q. Each magnetic 
Bloch band has one maximum, one minimum and two saddle points that have 
coinciding energies, and are linked by a separatrix. This is a consequence of 
Chambers relation [16]. For ip close to p/q, the q x q matrix valued function is 
quantized by imposing the commutation relation [£, 77] = 2ir i (<p — p/q). 

The semiclassical strategy leading to the qualitative description of the lip at 
ip = remains valid in this more general case, but the details are different. (For 
example, 7 does not take a universal value 1/2 but rather becomes a function 
of E and <p.) However, one still expects that the envelope of the magnetization 
is pinched linearly at the band edges and logarithmically at the separatrix. 
This gives all L b (p,,p/ q, j), j = 1, . . . ,q, a lip-like shape, a distorted version of 
Fig. 2. Getting explicit expressions for L b (/j,p/q,j) is a hard problem which lies 
beyond the scope of this paper. However, numerical support for our claims is 
given in Fig. 4. Here the magnetization as a function of the chemical potential 
is shown for the flux tp = 1 1 which is close to A. The three bands of ip = A 

•3+4(5 J 3 

split and give rise to rapidly oscillating magnetization with distorted lip-shaped 
envelopes. 

Band splitting occurs on all energy scales of the Hofstadter butterfly, because 
of its fractal nature. We have shown that band splitting is always accompanied 
by de Haas-van Alphcn oscillations in the magnetization. The fractality of 
the spectrum is reflected by the fact that the magnetic oscillations occur on 
arbitrarily small energy scale. The fingerprints of fractality are most evident 
for flux ratios tp which are irrational numbers that are well approximated by 
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rationals (e.g. Liouvillc numbers). For such ip the graph of the magnetization 
as a function of n, such as those shown in Figs. 2 and 4, is itself fractal, with 
magnetic oscillations bounded by lip- like envelopes appearing on infinitely many 
scales. 
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